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PARAMETRIC  STUDY  OF  THE  DYNAMIC  STABILITY  OF  TOWED  VESSELS 

Several  accidents  in  towing  operations  of  barges  or 
disabled  ships  in  restricted  and  open  waters  have  made 
necessary  the  investigation  of  the  course  keeping  stability 
of  towed  vessels.  In  this  work  a  non-linear  mathematical 
model  is  used  to  simulate  the  slow  surge,  sway,  and  yaw 
motions  of  a  vessel  towed  by  a  heavy  catenary  towline.  The 
effect  of  geometric  parameters  of  the  system  on  the 
stability  of  equilibrium  configurations  is  analyzed. 

It  is  shown  that  for  certain  choices  of  towing  system 
parameters,  dynamic  loss  of  stability  may  occur  which 
results  in  qualitatively  different  asymptotic  response.  The 
results  of  this  study  identify  regions  in  the  parameter 
space  that  lead  to  either  safe  operations  or  hazardous 
system  response. 
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I.  INTRODUCTION 

A.   BACKGROUND 

A  long  history  of  towing  accidents  resulting  in  loss  of 
life,  damage  to  property,  and  pollution  of  the  environment 
have  prompted  many  studies  into  the  dynamics  of  towing 
operations.  Of  primary  concern  was  the  motions  of  the  towed 
vessel  in  the  horizontal  plane  (yaw,  sway,  and  surge) . 
Excessive  and  unstable  motions  could  lead  to  collisions  and 
capsizing.  The  ability  to  predict  the  motion  of  a 
particular  towing  system  would  be  of  particular  benefit  to 
ship  designers  and  towing  operators,  by  identifying  those 
situations  where  the  towing  operations  would  be  the  safest, 
or  those  which  must  be  avoided. 

Previous  studies  at  the  University  of  Michigan  and 
elsewhere  had  developed  mathematical  models  and  numerical 
techniques  to  analyze  towing  dynamics,  and  had  identified 
those  parameters  which  are  of  primary  importance  to  the 
stability  of  the  towing  system.  The  linear  model  usually 
used  to  describe  ship  motions  [Ref.  1,  Chapter  7]  is 
inadequate  for  the  towing  problem.  Non-linear  models,  as  in 
[Ref.  2],  must  be  utilized  to  accurately  describe  the 
towing  system.  These  studies  had  identified  the  position  of 


the  towline  attachment  point  on  the  towed  vessel  and  the 
towline  tension  as  the  most  significant  (and  controllable) 
parameters  of  the  towing  system. 

B.   PROBLEM  CONDITIONS 

In  this  study,  computer  programs  developed  in  [Ref.  3] 
were  used  to  analyze  the  effect  of  different  parameter 
combinations  on  the  towed  stability  of  three  vessels.  These 
programs  use  bifurcation  to  identify  the  unstable  and 
stable  regions  of  the  parameter  space.  The  principal 
parameters  studied  were  (Fig.l): 

1.  longitudinal  position  of  the  towline  attachment  point 
forward  of  the  towed  vessel's  center  of  gravity,  xp; 

2.  athwartships  position  of  the  towline  attachment  point 
port  or  starboard  of  the  towed  vessel's  centerline, 

yp; 

3.  length  of  the  towline,  Lw. 

In  the  model  used  in  this  study,  unlike  [Ref.  2],  the 
towline  is  modeled  as  an  inextensible  caternary,  thus 
making  towline  tension  a  function  of  its  length.  The  model 
conditions  were: 

1.  speed  of  towing  vessel  of  2  knots; 

2.  towing  vessel  on  steady  course; 

3.  calm  seas,  no  wind;  i.e.,  no  external  environmental 
forces. 

Characteristics  of  the  towed  vessel  were  inputted  into  the 

programs  from  a  data  file  containing  hydrodynamic 

coefficients,  resistance  data,  towline  characteristics,  and 


TOWED  VESSEL 


TOWING  VESSEL 


Figure  1.   Problem  Geometry 


skeg  data,  if  applicable.  The  effect  of  asymmetrical  forces 
acting  on  the  towed  vessel,  such  as  the  presence  of  a 
propellor  or  an  environmental  force,  are  introduced  through 
a  bias  in  the  data  file.  All  dimensions  are 
nondimensionalized  with  respect  to  the  towed  ship's  length 
between  perpendiculars  (LBP) . 

Three  vessels  were  studied  (Fig. 2): 

1.  a  191  foot  barge  with  a  skeg,  with  no  propellor 
(i.e.,  no  bias) ; 

2.  a  1066  foot  tanker  with  no  skeg,  but  with  a  propellor 
(i.e.,  with  a  bias); 

3.  the  same  barge  as  in  1),  but  without  the  skeg  and 
with  a  propellor  (i.e.,  with  a  bias). 

Unlike  previous  studies,  this  work  includes  the  effect 
of  athwartship  position  of  the  towline  attachment  point  in 
the  stability  of  the  towing  system. 

Chapter  II  provides  background  into  the  problem 
formulation  and  stability  analysis  used  in  this  study. 
Chapter  III  presents  the  results  of  the  analysis  and 
discusses  some  practical  aspects  of  these  results.  Chapter 
IV  discusses  the  conclusions  which  can  be  made  from  the 
results  of  this  study  on  the  stability  of  the  towing  system 
and  the  use  of  the  techniques  used  herein. 
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*  Barge  with  skeg,  no  bias 
*  Barge  with  no  skeg,  with  bias 


Tanker  with  bias 


Figure  2 .   Body  Plans  of  Vessels  Studied 


II.  PROBLEM  FORMULATION  AND  METHOD  OF  APPROACH 

Slow  motions  of  a  towed  vessel  in  the  horizontal  plane 

are  described  by  a  system  of  six  nonlinear,  coupled, 

differential  equations.  [Ref.  4  and  5]  In  its  standard  form 

this  system  is 

1 
xl  =       ■~[F1(x1,X2,X3)  +  Tsurqe(x4,x5,x6) ] 
m  -  Xu 


1z    ~   Nr 

x2  ""         tF2  (xi'x2'x3^  +  Tsway  (x4 'X5'X6^ 
D 


Yr 
+     [F3 (x1,x2,x3)  +  xpTsway(x4,x5,x6)  ypTsurge(x4 ,x5,x6) ] , 


Nv 
x3  ~     tF2 (xl'x2'x3^  +  xpTsway^x4'x5'x6^ 


m  -  Y^ 
+        tF3 (xi'x2'x3)+xpTsway(x4'x5'x6^  ~ypTsurge (x4 ' x5 'x6^ 3 ' 


x4  =  x-^cosxg  -  x2sinx6  -  U, 


x5  =  x-^sinxg  +  x2cosx6, 


x6  ~  x3' 


where 

^surge  (x4  '  x5  '  x6^  =  ^x^x4  ,x5'x6^  cosx6+^y  ^x4  ,x5'x6^  sinx6' 

-r^sway  (x4  '  x5 '  x6^  =  ^x^x4  ,x5'x6^  sinx6~^y  ^x4 /X5'X6^  cosx6' 
D  denotes  the  known  quantity 

D  =  (m-Yv) (Iz-Nr)-YrNv, 

and 

2  3    i  2    i  2 

^1  (xi  /  x2  »  x3 )  ~xuxl+liXuuxl+1//6Xuuuxl+liXwx2  +  ^Xwux2xl 

2    i  2 

■  '2-"j-£.X3  + ^X— r uX3Xl     (*yr  ^J  X2X3~*~    rvuXlX2X3  ' 


2  3  2 

F2 (xx , x2 , x3 ) =YQ+Y0ux1+Youux1+Yvx2+l/6YvvvX2+^YvrrX2X3 

2  3   i       2 

+Yvuxlx2+^Yvuux2xl+  (^r~mxl^  x3+^-/^^rrrx3+'5^rwx3x2 


j        2 
+Yrux3xl"l"iYruux3xl ' 


F3 (xx , x2 , x3 ) =N0+N0ux1+Nouux1+Nvx2+l/6Nvvvx2+^Nvrrx2X3 

2  3  2 

+Nvux1x2+^NvuuX2x1+Nrx3+l/6Nrrrx3+^NrvvX3X2 

x        2 
+NruX3x1+^Nruux3x1. 

In  the  above  equations,  x-^  denotes  the  sway  velocity  in 
surge  (longitudinal  motion)  of  the  towed  vessel,  x2  the 
velocity  in  sway  (lateral  motion) ,  x3  the  angular  velocity 
in  yaw  (turning  motion  about  the  vertical  axis) ,  x4  and  x5 
the  coordinates  of  the  center  of  gravity  of  the  towed  vessel 
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with  respect  to  an  (x,y) -coordinate  system  moving  with  the 
towing  vessel,  and  x6  the  towed  vessel  yaw  angle.  Further,  U 
is  the  steady  towing  vessel  velocity  in  the  x-direction,  xp 
and  yp  are  the  coordinates  of  the  towline  connection  point 
on  the  towed  vessel  with  respect  to  an  (X,Y) -coordinate 
system  with  its  origin  at  the  towed  vessel  center  of 
gravity,  and  Rx,  Ry  are  the  towline  restoring  forces.  The 
towing  system  configuration  and  notation  conventions  are 
shown  in  Figure  1.  Expressions  for  F-^,  F2 ,  F3  are  derived  by 
Taylor  expansion  in  terms  of  the  relative  velocities  x^,  x2 , 
x3  of  the  towed  vessel  with  respect  to  the  water.  In 
nonlinear  analysis  terms  up  to  third  order  are  used  whereas 
terms  beyond  third  order  and  second-  and  higher-order 
acceleration  terms  are  usually  neglected.  Subscripts  u,  v,  r 
indicate  derivative  of  force-moment  component  with  respect 
to  Xn,  x2 ,  x3  respectively,  and  subscript  c  indicates 
propel lor  dependent  terms.,  which  represent  a  source  of 
system  asymmetry.  These  terms  are  zero  in  the  absence  of  a 
propellor.  Terms  Xabc,  Yaj->c,  Nabc,  where  a,  b,  c  are  dummy 
independent  variables  representing  u,  v,  r,  are  usually 
called  slow  motion  derivatives.  In  unsteady  reference 
motion,  slow-motion  derivatives  are  considered  as  functions 
of  the  frequency  of  motion.  In  our  study  of  slowly  varying 
reference  motions,  we  assume  that  slow  motion  derivatives 
are  time  independent.  This  is  a  good  approximation  for  ships 
with  usual  hull  shapes  and  moderate  speeds. 
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Rx  and  Ry  denote  restoring  forces  from  the  towline,  and 
for  a  quasistatic  towline  response  they  are  expressed  as 
implicit  functions  of  x4 ,  x5,  x6.  In  this  study,  the  model 
used  for  the  towline  is  that  of  an  inextensible  heavy 
caternary  with  nonlinear  force-displacement  characteristics 
as  given  in  [Ref.  3]. 

In  compact  notation  the  above  system  of  six  ordinary 
differential  equations  is  denoted  as 

x=f(x)  (1) 

where  x  and  f  are  six  dimensional  vectors.  To  analyze  the 
stability  properties  of  (1) ,  the  first  step  is  to  identify 
the  equilibrium  configuration  of  the  system.  For  this  we 
have  to  solve  a  system  of  six  nonlinear,  coupled  algebraic 
equations 

f(x)=0  (2) 

where  x  denotes  an  equilibrium  configuration.  It  can  be 
shown  [Ref.  5]  that  system  (2)  has  at  most  three  solutions 
in  x  corresponding  to  three  distinct  equilibrium  positions. 
In  this  study  we  concentrated  our  efforts  on  one  of  these 
equilibrium  positions,  namely  the  one  which,  in  the  absence 
of  a  bias  in  the  system,  corresponds  to  the  towed  vessel 
being  located  directly  astern  of  the  tow-tug.  This  is  the 
most  interesting  in  applications.  Having  computed  x,  its 
stability  properties  can  be  established  as  follows: 
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Linearization  of  (1)  around  x  leads  to  the  linear  system 

z=Az  (3) 

where  z  represents  the  excursion  from  the  equilibrium  x,  and 
A  is  a  constant  6x6  matrix.  If  all  eigenvalues  of  A  have 
negative  real  parts,  then  x  is  stable;  otherwise  it  is 
unstable. 

In  this  study,  we  performed  parametric  analysis  of  the 
central  equilibrium  in  therm  of  towline  length  Lw,  and  the 
towing  point  coordinates  with  respect  to  the  center  of 
gravity  of  the  towed  vessel,  xp  and  yp.  These  parameters  can 
be  easily  changed  before  or  during  towing  operations  and  can 
provide  a  means  of  passive  control  of  the  towing  system. 
Parameter  Lw  directly  affects  the  amount  of  tension 
developed  by  the  towine.  Parameter  xp  is  directly  related  to 
the  towline  restoring  force  and  moment.  A  small  value  for  xp 
may  not  be  able  to  provide  adequate  restoring  moment  and  may 
not  guarentee  system  stability.  On  the  other  hand  a  very 
large  value  for  xp  may  result  in  over-compensation  and 
therefore  instability.  Nonzero  yp  values  result  in  a  source 
of  asymmetry  introduced  in  the  system.  For  a  biased  system 
(for  example  due  to  the  presence  of  a  propellor  or 
environmental  forces) ,  it  should  be  expected  that  an  extra 
appropriate  bias  introduced  via  a  nonzero  yp  helps 
counteract  the  effect  of  the  former  bias,  and  hence  improve 
stability. 
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The  particular  equilibrium  position  will  lose  its 
stability  when  an  eigenvalue  of  the  A  matrix  in  (3)  changes 
its  sign  from  negative  real  part  to  positive  real  part.  The 
case  when  a  real  eigenvalue  crosses  zero  has  ben  analyzed  in 
detail  in  [Ref.  5],  This  corresponds  to  a  static  loss  of 
stability  with  generation  of  additional  equilibrium 
positions  in  the  form  of  solution  branching.  In  this  study 
we  focussed  our  attention  on  the  case  when  a  complex 
conjugate  pair  crosses  the  imaginary  axis.  This  corresponds 
to  a  Hopf  bifurcation:  the  particular  equilibrium 
experiences  a  dynamic  loss  of  stability  and  the  system 
begins  to  oscillate.  The  resulting  periodic  solutions  can  be 
stable  or  unstable,  but  at  any  rate,  such  a  situation  is 
hazardous  and  should  be  avoided  during  towing  operations. 
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III.  RESULTS  AND  DISCUSSION 

A.   BARGE  WITH  SKEG 

The  first  vessel  to  be  studied  was  an  unpowered  barge 
with  a  skeg  aft.   Since  there  is  no  propellor  to  introduce 
a  bias,  the  barge  has  athwartship  symmetry. 

1.   Figure  3:  Critical  Real  Part  vs.  xp 

Program  T0WBIF1  calculates  eigenvalues  for  specific 
Lw  and  yp,  creates  a  file  for  each  of  six  real  and  six 
imaginary  parts,  and  creates  a  separate  file  containing  the 
largest  real  part.   The  real  part  with  the  largest  value  is 
the  critical  indicator  of  the  system's  stability:  if  it  is 
greater  than  zero,  the  system  will  be  unstable;  if  less 
than  zero,  the  system  will  be  stable. 

Figure  3  shows  plots  for  the  critical  parts  for 
yp=0.05  and  three  values  for  Lw.   The  region  where  the  plot 
is  greater  than  zero  indicate  that  range  of  xp  where  the 
system  is  unstable.  For  example,  for  Lw=0.5  the  critical 
real  part  is  greater  than  zero  for  the  range  of  xp=0.18  to 
xp=0.48,  so  the  system  is  unstable  within  this  range. 

Note  that  as  Lw  increases,  the  unstable  range 
becomes  smaller,  until  for  Lw=3.0  there  is  no  region 
greater  than  zero.   Therefore,  the  system  will  be  stable 
for  all  values  of  xp;  i.e,  the  barge  should  exhibit  no 
unstable  motion. 
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2.  Figure  4:  yp  vs.  xp,  Lw  as  Parameter 

Program  T0WBIF2  does  the  same  calculations  as 
T0WBIF1  over  a  range  of  values  of  yp  with  a  given  Lw, 
instead  of  a  single  value  of  yp  and  Lw.   In  essence,  Figure 
3  represents  a  cut  of  Figure  4  at  a  single  value  of  yp  and 
Lw.  Unlike  T0WBIF1,  T0WBIF2  writes  a  point  only  where  the 
critical  real  part  changes  sign.   When  plotted,  these 
produce  a  curve  delineating  stable  and  unstable  regions  of 
the  parameter  space.   This  is  the  point  of  the  process;  we 
are  more  interested  in  finding  what  parameters  produce 
stable  or  unstable  system  than  the  actual  results  of  the 
equations  of  motion. 

Recalling  Figure  3,  The  area  inside  the  curves 
represent  the  unstable  region.   For  example,  for  Lw=0.5, 
the  system  is  unstable  for  all  values  of  yp  within  the 
range  xp=0.2  to  xp=0.5.   Increasing  Lw  first  decreases  the 
unstable  range  of  xp  for  high  yp,  then  decreases  the 
unstable  range  of  yp.  For  large  Lw  (>4.0),  the  unstable 
range  virtually  disappears. 

3.  Figure  5:  Lw  vs.  xp,  yp  as  parameter 

Program  T0WBIF3  performs  the  same  calculations  as 
T0WBIF2 ,  but  with  Lw  as  the  ordinate  and  yp  as  the 
parameter.   Thus  Figure  5  provides  the  same  information  as 
Figure  4  but  with  a  different  perspective. 
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As  in  Figure  4,  the  unstable  region  is  inside  the 
curves.   It  clearly  shows  how  increasing  Lw  decreases 
the  unstable  range  for  a  constant  yp,  as  was  evident  in 
Figure  4.   It  also  shows  that,  for  constant  Lw  greater  than 
about  0.7,  increasing  yp  also  decreases  the  unstable 
region.   In  the  narrow  range  of  Lw  from  0.2  to  0.7, 
increasing  yp  increases  the  extent  of  the  unstable  range  of 
xp.  This  effect  is  apparent  in  Figure  4,  but  more 
dramatically  presented  in  Figure  5.   It  would  appear  that 
using  two  views  of  the  data  would  emphasize  aspects  of  the 
curve  that  may  be  overlooked  with  one  view. 

Since  positive  values  of  yp  represent  port  side 
placement  of  the  towline  attachment  point,  and  negative 
values  starboard  side  placement,  both  positive  and  negative 
values  for  yp  were  studied.  As  expected  from  the  port- 
starboard  symmetry  of  the  barge,  curves  for  positive  and 
negative  values  of  yp  were  identical,  and  only  positive 
values  were  presented  here. 

From  an  operational  point  of  view,  one  may  conclude 
from  these  curves  that  for  the  unpowered  barge,  placing  the 
towline  on  an  attachment  point  to  either  side,  as  far 
forward  as  possible,  will  make  the  tow  stable  for  the 
greatest  range  of  towline  length,  but  the  towline  should  be 
kept  no  shorter  than  the  length  of  the  barge. 
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B .   TANKER 

The  second  vessel  studied  was  a  tanker  typical  of  those 
now  in  service.  The  effect  of  the  tanker's  propellor  makes 
the  hull  asymmetrical;  this  effect  is  represented  by  a  bias 
included  in  the  tanker  data  file. 

1.  Figure  6:  Critical  Real  part  vs  xp 

Figure  6  plots  data  generated  from  T0WBIF1  with 
yp=0.10  and  two  values  for  Lw.  These  plots  show  the  stable 
region  to  be  between  the  two  zero  values  for  the  curve. 
Note  that  the  stable  region  becomes  smaller  with  increasing 
Lw. 

Note  also  that  both  curves  are  discontinuous  in  their 
slopes.  The  critical  real  parts  file  is  a  composite  of 
several  results  files,  each  of  which  is  critical  over  a 
certain  range.  Each  results  file  forms  a  smooth  curve; 
thus  the  curves  plotted  on  each  T0WBIF1  figure  may  be 
combinations  of  the  critical  section  of  several  results 
files. 

Finally,  note  that  the  stable  region  occurs  over  a 
narrow  range  of  xp,  unlike  the  barge  with  skeg  discussed 
earlier. 

2.  Figure  7:  Lw  vs.  xp,  yp  as  Parameter 

Figure  7  was  produced  from  data  generated  by  T0WBIF3 
for  positive  values  of  yp.  As  was  shown  in  Figure  6,  the 
stable  region  is  inside  the  curves.   The  vertical  line  at 
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xp=0.404  is  a  common  crossing  point  for  all  curves.  Each 
curve  is  formed  by  two  cusps,  with  the  upper  cusp 
dominating  with  decreasing  yp.   Each  cusp  is  the  plot  of 
different  critical  pair  of  eigenvalues;  the  "nose"  in  the 
curves  is  the  point  where  they  intersect. 

Note  how  the  stable  region  gets  smaller  with 
decreasing  yp  for  Lw  less  than  0.7,  for  example,  with 
yp=0.0  and  Lw  greater  than  1.0,  there  is  a  very  narrow 
range  of  xp  where  stability  can  be  assured. 

3.  Figure  8:  Critical  Real  Part  vs.  xp 

T0WBIF1  was  again  used  to  form  Figure  8,  this  time 
with  one  value  for  Lw  (Lw=0.6)  and  three  negative  values 
for  yp.  The  negative  yp  curves  pass  from  stable  to  unstable 
regions,  with  the  stable  ranges  for  xp  getting  smaller  as 
xp  becomes  more  negative. 

As  in  Figure  6,  the  curves  are  composites  of  those 
results  curves  which  are  critical  over  a  particular  range 
of  xp. 

4.  Figure  9:  Lw  vs.  xp,  Negative  Values  of  yp  as 
Parameters. 

Figure  9  data  was  generated  from  T0WBIF3,  with  yp=0.0 
curve  included  to  provide  continuity  with  Figure  7. 

The  stable  region  gets  smaller  as  yp  decreases  from 
0.0.  At  yp=-0.10  the  "nose"  between  upper  and  lower  cusps 
appears  to  be  tipping  up,  with  the  region  inside  the  "nose" 
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being  stable.   The  curves  plotted  in  Figure  8  were  formed 
using  a  value  of  Lw  which  cut  through  this  nose,  thus 
forming  the  sinuous  curves  which  pass  in  and  out  of  the 
stable  region.   Note  that  for  the  most  negative  values  of 
yp,  the  cusps  have  disappeared,  and  the  stable  range  of  xp 
is  slightly  increasing. 

5.   Figure  10:  Lw  vs  xp,  Negative  Values  of  yp  as 
Parameters. 

Figure  10  is  a  "close-up"  of  Figure  9,  focusing  on 
what  is  happening  around  yp=-0.10.  The  lower  cusp  tips  up 
and  merges  with  the  upper  to  form  a  single  curve.  Note  how 
the  stable  range  of  xp  virtually  disappears  for  Lw  greater 
than  0.5  for  yp=-0.10  and  -0.11.  As  was  seen  in  Figure  9, 
stable  range  for  xp  for  Lw  greater  than  0.5  reappears  with 
yp  less  than  -0.15. 

C.   BARGE  WITHOUT  SKEG 

The  third  vessel  was  a  self-propelled  version  of  the 
barge  studied  in  Section  A  (not  under  power  during  tow) , 
but  without  the  skeg.  As  with  the  tanker,  the  presence  of 
the  propellor,  simulated  by  a  bias  in  the  data  file, 
introduces  port-starboard  asymmetry. 

1.   Figure  11:  Critical  Real  Part  vs  xp  for  Lw=1.5 
Figure  11  shows  curves  for  three  values  of  yp 
(greater  than  zero,  zero,  and  less  than  zero)  and  one  value 
of  Lw.  The  curves  show  the  stable  range  of  xp  to  be  between 
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the  zero  crossing  points  of  the  critical  real  parts,  as  in 
the  tanker  case.  Also  similarly  to  the  tanker,  the  stable 
region  increases  with  increasing  yp.  These  results  are 
opposite  to  the  propellor-less  barge  with  the  skeg. 

2 .   Figure  12 :  Lw  vs  xp 

Figure  12  dramatically  shows  how  decreasing  yp 
reduces  the  stable  region.  The  vertical  line  at  xp=0.14  was 
common  to  all  values  of  yp  greater  than  and  equal  to  zero. 
For  values  of  yp  less  than  zero,  the  smooth  shape  of  the 
curve  is  apparent. 

The  tanker  and  self-propelled  barge  cases 
dramatically  demonstrate  the  effect  that  a  bias,  like  a 
propellor,  can  introduce  to  the  stability  of  the  system. 

D.   PRACTICAL  OBSERVATIONS 

Analysis  of  the  graphs  suggests  some  general  principles 
which  may  be  applied  when  conducting  slow  speed  towing 
operations  with  the  vessels  discussed  in  this  chapter. 
While  these  principles  are  of  course  not  generally 
applicable  to  all  vessels,  they  illustrate  how  the  analysis 
techniques  employed  in  this  work  can  be  applied  to  other 
vessels. 

For  the  unpowered,  symmetric  barge  with  a  skeg,  the 
operator  should  have  the  towline  attachment  point  as  far 
out  to  either  side  as  possible  and  the  towline  as  long  as 
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practical.   The  attachment  point  can  then  be  placed  at  any 
location  forward  of  the  center  of  gravity  with  stability 
assured.   Conversely,  if  the  attachment  point  must  be  on 
the  centerline,  placing  it  as  far  forward  as  possible 
(about  half  the  barge's  length  forward  of  the  center  of 
gravity)  will  assure  stability  for  all  towline  lengths. 
For  vessels  with  an  asymmetrical  bias  (e.g.,  with  a 
propellor) ,  but  without  skegs,  the  attachment  point  needs 
to  be  as  far  to  the  biased  side  as  possible  (in  the  cases 
of  the  tanker  and  self-propelled  barge,  the  +yp  or  port 
side)  and  placed  forward  of  the  center  of  gravity  the 
distance  indicated  on  the  graph  for  all  towline  lengths. 
Placement  of  the  attachment  point  on  the  opposite  side  (in 
the  cases  studied,  the  starboard  side)  will  virtually 
assure  the  system  to  be  unstable. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

This  study  highlighted  the  effect  of  athwartship 
position  of  the  towline  attachment  point.   The  common 
assumption  among  ship  operators  prior  to  this  reseach  held 
that  placing  the  towline  on  the  centerline  on  the  foremost 
point  of  the  towed  vessel  would  create  the  most 
advantageous  towing  situation.   Studies  such  as  [Ref.  2] 
have  shown  that  towing  stability  can  be  dependent  on  the 
longitudinal  placement  of  the  attachment  point.   This 
research  has  shown  that  for  certain  conditions,  attaching 
the  towline  off  the  centerline  can  also  improve  towed 
stability.   The  optimum  towing  configuration  requires  a 
combination  of  all  three  parameters  -  longitudinal  and 
athwartship  placement  of  the  towline  attachment  point,  and 
towline  length. 

The  bifurcation  technique  used  in  this  study  can  be 
used  to  produce  stability  information  useful  to  ship 
designers  and  towing  operators.  Stability  information  can 
be  assembled  into  a  convenient  graphical  form  that  clearly 
defines  the  regions  of  stable  and  unstable  operation  based 
on  the  parameters  the  operator  has  the  most  control  over  - 
the  placement  and  length  of  the  towline. 
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For  the  ship  designer,  this  technique  can  be  useful  in 
determining  the  implications  particular  design  decisions 
would  have  on  the  vessel's  performance  under  tow. 
Depending  on  the  vessel's  use,  adjustments  to  the  design 
can  be  made  to  improve  towing  stability,  or  the  customer 
can  be  forewarned  to  avoid  certain  kinds  of  operations. 
Since  nearly  all  vessels  are  towed  at  some  time,  towing 
performance  should  be  analyzed  for  all  vessels. 

For  the  towing  operator,  this  technique  can  provide 
readily  available  information  about  how  a  particular  vessel 
will  respond  under  tow.   The  operator  can  then  adjust  the 
towing  parameters  (i.e.,  placement  the  attachment  point 
and/or  length  of  the  towline)  so  the  tow  will  be  in  its 
most  stable  condition,  or,  if  unavoidable,  know  that  a 
particular  towing  situation  will  be  potentially  dangerous 
and  make  preparations  to  deal  with  it. 

Since  ship  data  is  inputted  through  a  data  file,  the 
towed  performance  of  any  vessel  can  be  analyzed  with  this 
method,  including  structures  such  as  offshore  oil 
platforms.   Existing  vessels  can  be  analyzed,  as  well  as 
different  loading  conditions. 
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Two  principal  disadvantages  are  associated  with  this 
technique: 

1.  The  programs  are  dependent  on  the  quality  of  the 
data  provided.  Determining  hydrodynamic  coefficients 
and  resistance  data  requires  tow  tank  experiments  and 
analysis,  and  are  not  obtained  for  most  vessels; 

2 .  The  programs  require  large  amounts  of  computer  time 
and  memory  to  run,  which  may  not  be  available  or  too 
costly  for  potential  users,  especially  to  run 
extensive  "what-if"  scenarios.  This  problem  may  be 
alleviated  as  more  inexpensive,  high  speed,  high 
capacity  micro-  and  personal  computers  become 
available. 


B.  RECOMMENDATIONS 

This  study  was  done  for  only  one  set  of  conditions. 
Further  research  can  be  done  in  determining  the  effect  of 
varying  conditions,  such  as  different  speeds  or  maneuvering 
by  the  towing  vessel,  on  the  stability  of  the  tow.  External 
forces  are  modelled  by  the  bias  in  the  data  file.  A 
systematized  method  of  introducing  biases  into  the  data 
would  enable  the  analysis  of  the  effect  of  environmental 
conditions  on  the  towing  system. 

Further  work  should  be  conducted  to  improve  the  "user- 
friendliness"  of  the  programs.  As  currently  configured,  the 
programs  must  be  run  instructively,  and  graphics  produced 
offline.  This  is  a  time  consuming  process  which  does  not 
use  the  full  capabilities  of  either  the  programs  or  the 
graphics  capabilities  of  the  mainframe.  Program 


32 


improvements  should  focus  on  streamlining  computations  and 
user  interaction,  and  incorporating  graphics,  with  the  goal 
of  making  it  available  as  a  ship  design  tool. 
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APPENDIX 


Driver  programs  used  in  this  thesis  are  shown  here 
Subroutines  can  be  obtained  by  contacting: 

Prof.  F.A.  Papoulias 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 

Monterey,  CA  9394  3 

FILE:  T0WBIF1   FORTRAN   Al 


PROGRAM  T0WBI1 


BIFURCATION  ANALYSIS  OF  TOWING  SYSTEMS 
PARAMETER  DEPENDS  ON  IPAR 


IPAR  =  1 


XP 
YP 
LW 


IT  NEEDS  SUBROUTINES  FROM  TOWING. FTN 

IMPLICIT  DOUELE  PRECISION  (A-H.O-Z) 

DOUELE  PRECISION  MA33P .NVD.NV.NRD. IZZ.NR.LB.LEN.LW, 

1  NO . NOU . NOUU • NVVV . NVRR . NVDD . NVU . NVUU . NRRR . NR W . 

2  NRDD . NRU . NRUU . ND • NDDD . NDVV . NDRR . NDU . NDUU . NVRD 

DIMENSION  IVK6).AC6.6).VV(3).X(6).WR((>).WI(6).ZC6,6).SV2C6) 

COMMON/INTGR/ISKEG.NREDP . ITYS . ID. IFDS. I  STAB. IPROP 

COMMON/SPAR /MASSP . LW . XPP . YPP . LB 

CCMMON/SURGE/SU( 7 ) 

COMMON/XSURG/XU.XUU.XUUU 

CCMMON/SWAY/SW( 15 ) 

COMMCN/YAW/YA( 16 ) 

C0MM0N/MTER/VCAR.RH0.ABS-C0N1 .CON2 

COMMON/RES  I ST/VELKO).  RES  I  (40) 

COMMON/VELE/UEL( 100) 

COMMCN/POSTN/X1 .Yl .21 

COMMCN/GEOM/AL . R W . G . AET . HW . HW 1 

COMMON/PROP/ALE. P. EY.DI A. ANIU 

COMMCN/CTNR/XC ( 99 ) . YC ( 99 ) . ZC( 99 ) . TC<  99 ) 

COMMON/INT1/IC 

COMMON/DOC/UC. ALPHA 

COMMON/UEPT/RLX . RLY .RLZ 

COMMON/SLOPE /PDRXX-PDRXY.PDRYX.PDRYY 

COMMON/SLAN/RXX6.RYY6  .RXX.RYY 

OPEN  (UNIT  =  35.FILE  =  ,BARGE:i  ,STATUS= 'OLD'  ) 
CFE'J  (UNIT  =  1  ,FILE  =  'RE30'  .STATUS=  "NEW  ) 

OF  EN  (UNITrll.FILE^RESlR'  ,STATU3='NEW  ) 

OPEN  (UNIT=12.FILE='RES2R- . STATUS = • NEW ) 

OPEN  CUNIT=13.FILE='RESIR' .STATUS' 'NEW ) 

OPEN  CUNIT=1«.FILE='RES4R' .STATUS^ ' NEW  ) 

OFEN  (UNITrlB-FlLE^RESBR1  ,STATUS='NEW  ) 

OPEN  (UNIT^e.FILEi'RESbR'  ,STATUS='NEW  ) 


TOW00010 
TOW00020 
TOW00030 
TOW00040 
TOW00050 
TOW00060 
TOW00070 
TOW00C8O 
TOW00090 
TOW00100 
TOW00110 
TOW00120 

towooi:o 

TOW00K0 
TOW00150 
TOW00160 
TOW00170 
TOW00180 
TOW00190 
TOW0C200 
TOW00210 
TOW00220 
TOW00230 

towoo:<;o 

TOW00250 
TOW00260 
TOW00270 
TOW00280 
TOW00290 
TOW00300 
TOW00310 
TOW00320 
TOW00330 
TOW003«0 
TOW00350 
TOW00360 
TOW00370 
TOW00380 
TOW00390 

towoosoo 
towoosio 

TOW00420 

towoo<;30 

TOW00440 
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OPEN  (UNIT=21 .FILE='REC1I 

OPEN  (UNIT=22-FILE='RES2I 

OPEN  (UNIT=23.F1LE='RES3I 

OPEN  (UNIT=24.FILE='RES4I 

OPEN  (UNIT=25.F1LE='RE35I 

OPEN  (UNIT=26-FILE='RES6I 

CALL  INPUTD(IO) 
VCAR  =  VCARM.689D0 
MATT  =0 
IFLOW=l 


WRITE  (».1001) 

READ   («■»)     IPAR 

WRITE  C  »  .  1 002  ) 

READ   («.«)     A1.A2 

WRITE  (".1003) 

READ   (",«)     NUM1 

WRITE  («.1005) 

READ   (».»)     1KB 

WRITE  t  »  .  1 006  ) 

READ   (».«)     NEOL 

IF  (IKB.GT.NECL)  GO  TO  500 

DO  1  I  =  1 . MUM  1 

WRITE  (».2001)  I.NUM1 
AA=A1'(A2-A1 )»( 1-1 J/INUM1-1 ) 
IF  (IPAR.EO.l)  XPP=AA 

IF  (IPAR.EO.2)  YPP=AA 
IF  (IPAR.EQ.3)  LW  =AA 

AL  =LW-LB»0.30«;8D0 

ALE=AL 

CALL  STABIL(IW.W.ISOL) 


•  ,STATUS='NEW'  ) 

•  .STATUS^'NEW  ) 

•  ,STATUS='NEW  ) 
'  .STATU3='NEW'  ) 
' ,STATU3='NEW' ) 
• .STATUS='NEW' ) 


SET  V=VV(K)  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 

IF  (IW.NE.NEOL)  GO  TO  1 
V  =  VV( 1KB) 

if  (dabs(v).gt. 1 .do)  stop  1111 

call  eouilb(v.x.res-rx.ry) 

call  linear(x.res.a.rx.ry) 

call  rg(6 -6.a.wr.wi.matz.z. iv1.sv2. ier1 ) 

if  (ieri .ne.o)  stop  22:2 

call  degstb(deos.wr) 

WRITE  (l.)0)  AA.DEOS 
DO  11  J=l .6 
JR=10*J 

WRITE  (JR. 10)  AA.WR(J) 
JI=20*J 

WRITE  (JI  ,  10)  AA.WKJ) 
11    CONTINUE 
1  CONTINUE 
500  STOP 
10  FORMAT  (2D20. 10) 


1001  FORMAT  C 
1 


1002  FORMAT  ( 

1003  FORMAT  ( 

1005  FORMAT  ( 

1006  FORMAT  ( 
2001  FORMAT  (215) 

END 


ENTER  1  :  XP  VARIATION'./, 

2  :  YP  VARIATION' ,/. 

3  :  LW  VARIATION' ) 
ENTER  PARAMETER  RANGE') 
ENTER  NUMEER  OF  INCREMENTS') 
ENTER  EQUILIBRIUM  NUMBER') 
ENTER  ESTIMATED  NO.  OF  EQUILIBRIA') 


TCW00<;50 
TCW00460 
TOW00<W0 
TOW00480 
TCW00<1°0 
TOW00500 
TOW00510 
TOW00520 
TOW00530 
TCW00540 
TOW00550 
TOW00560 
TOW00570 
TCW00580 
TCW00590 
TOW00600 
TCW00610 
TOW00620 
TOW0O630 
TCW00640 
TOW00650 
TCW00660 
TOW00670 
TCU00680 
TOW00690 
TCW00700 
TOW00710 
TCW00720 

TOW00730 
TOW007G0 
TCW00750 
TCW00760 
TOW00770 
TCW00780 
TCW00790 
TOW008C0 
TCW00810 
TOW00820 
TCW00830 
TOW008<;0 
TCW00850 
TOW00860 
TOW00870 
TOW00880 
TOW00890 
TOW00900 
TOW00910 
TOW00920 
TOW00930 
TOW009<;0 
TOW00950 
TCU00960 
TOW00970 
TCW00980 
TCW00990 
TCW01000 
TOW01010 
TOW01020 
TOW01030 
TOW01060 
TOW01050 
TOW01060 
TOW01070 
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FILE:  T0WBIF2   FORTRAN   Al 


PROGRAM  TOWEI2 
PROGRAM  TOWBIF.FTN 

BIFURCATION  ANALYSIS  OF  TOWING  SYSTEMS 

PARAMETERS  ARE:  Xp.  Yp 

IT  NEEDS  SUBROUTINES  FROM  TOWING. FTN 

USER  DEPENDENT  SUBROUTINES: 

DEGSTB  =  CURVES  ENCLOSING  REGION  II  OF  FIGURE  13 

(SUBROUTINE  DEGSTB  IS  IN  TOWING. FTN) 
DS1     =  CURVES  ENCLOSING  REGION  V  OF  FIGURE  13 

(SUBROUTINE  DS1  IS  IN  SPMBIF) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  MASSP .NVD.NV .NRD. IZZ .NR ,LB .LEN.LW, 

1  NO-NCU.NOUU.NWV.NVRR.NVDD.NVU.NVUU.NRRR  .NRVV, 

2  NRDD.NRU.NRUU.ND.NDDD.NDVV.NDRR .NDU .NDUU.NVRD 

DIMENSION  IV1(6).A(6.6).VV(3).XU).WR(6).WI(6).Z(6.6).SV2(6) 

COMMON/ 1NTGR/ 1 SKEG.NREDP . ITYS. ID. IFD3. 1  STAB. I  PROP 

COMMON/SPAR/MASSP . LW . XPP . YPP . LB 

COMMCN/SURGE/SU( 7 ) 

COMMCN/XSURG/XU.XUU.XUUU 

CCMMCN/3WAV/SW( 15 ) 

COMMON/YAW/YAt 16) 

COMMCN/MTER/VCAR.RHO.ABS.CON1 ,CON2 

COMMON/R  ES I  ST/ VEL (40). RES  1(40) 

COMMON/VELE/UEL (100) 

CCMMCN/POSTN/X1 .Yl ,21 

COMMON/GEOM/ AL . RW . G . AET . HW . HW1 

COMMON/PROP/ALE. P. EY.DI A. AN  I U 

COMMCS/CTNR/XC(99).YC(99).ZC(99).TC<99) 

COMMON/INT1/IC 

COMMGN/DOC/UC • ALPHA 

COMMON/UEPT/RLX. RLY . RLZ 

COMMCN/SLOPE/PDRXX.PDRXY.PDRYX.PDRYY 

COMMON/SLAN/RXX6 .RYY6 .RXX.RYY 

OPEN  (UN^iB-FILEi-BSKEGZ1 ,STATUS= 'OLD' ) 
OPEN  (UNIT=11 -FILE^'RESIR1 . STATU3= ' NEW ) 
OPEN  (UNIT=12.FILE='RES2R' . STATU3= ' NEW' ) 
OPEN  (UNIT=13.FILE='RES3R' . STATUS =' NEW ) 
OPEN  (UNIT=14.FILE=,RES'R- ,STATU3=,NEW ) 
OPEN  (UNITilB.ciLEs'RESBR1  .3TATUS='NEW  ) 
OPEN    (UNIT^lS.FILE^'RESbR' ,3TATUS='NEW ) 

CALL  INPUTD(IO) 
VCAR  =VCAR«1 .689D0 
AL    =LW«LB"0.3048D0 
ALE   =AL 
MATZ  =0 
IFLOW=l 
EPS   =l.D-5 
ILMAX=1500 


WRITE  (' 

READ  (' 

WRITE  (' 

READ  (' 

WRITE  (' 

READ  (' 

WRITE  (' 

READ  !' 


1001  ) 

«)     A1.A2 

1002) 

«)     NUM1 

1003) 

»)     B1.B2 

1004) 

«)     NUM2 


TOW00010 
TOW0002C 
TOW0003C 
TOW00040 
TOW000S0 
TOWOOOcO 
TOW00070 
TOW00080 
TOW00C90 
TOW00100 
TOW00110 
TOW00120 
TOW00130 
TOW00140 
TOW00150 
TOW00160 
TOW00170 
TOW00180 
TOW00190 
TOW00200 
TOW00210 
TOW00220 
TOW00230 
TCW0024O 
TOW00250 
TOW00260 
TOW0C270 
TOW00280 
TOW00290 
TOW00300 
TOW00310 
TOW00320 
TOW00330 
TOW00340 
TOW00350 
TOW00360 
TOW00370 
TOW00380 
TOW00390 
TOW00400 
TOW00M0 
TCW00420 
TOW00430 
TOW004  40 
TOW00450 
TOW00460 
TOW00470 
TOW00480 
TOW00490 
TOW00500 
TOW00510 
TOW00520 
TOW005S0 
TOW00540 
TOW00S50 
TOW00560. 
TCW00570 
TOW00580 
TOW00590 
TOW00600 
TOW00610 
TOW0C620 
TCW00630 
TOW00t>40 
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WRITE  C.10C5) 
READ   (».«)     1KB 
WRITE  C.1006) 
READ   («.«)     IDG 
DO  1  I=1.NUM1 

WRITE  <«.2001)  I.NUM1 

YPP=A1»(A2-A1 )«(I-1 )/(NUMl-l ) 


TOW00650 
TOW00660 
TOW00670 
TOWOO08O 
TOW00690 
TOW00700 
TOW00710 
TOW00720 


XPP=B1 

CALL  STABILC  IW.W.ISOL) 

SET  V=VV(K)  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 

V=WC!KB) 

IF  (DABS(V).GT.l.DO)  STOP  1111 

CALL  EOUILB(V.X.RES.RX.RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RG(6,6.A.WR.WI .MATZ.Z.IV1 -SV2.IER1) 

IF  (1ER1 .NE.O)  STOP  2222 

IF  (IDS.EO.l)  CALL  DEGSTBtDECS.WR) 

IF    (ID3.E0.2)    CALL    D3KDE0S-WR) 

DE0300=DE03 

XPOO  =XPP 

L     =0 

DO  2  J=2.NUM2 

WRITE  («,«)  J 

XPP=B1*(B2-B1 )«( J-l 1/CNUM2-1 ) 

CALL  STAEILCIVV.VV.ISOL) 

SET  V=VV(K)  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EOUILIERIUM 

V=VV( 1KB) 

IF  (DAB3(V).GT.1.D0)  STOP  1111 

CALL  EOUILB(V.X.RES.RX.RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RG(6.6.A.WR.WI .MATZ.Z.IV1 .SV2.IER1 ) 

IF  (IER1 .NE.O)  STO=  2222 

IF  (ID3.EQ.1)  CALL  DEG3TB (DEOS . WR ) 

IF    (IBS.EQ.2)    CALL   DSKDEOS.WR) 

DE0SNN=DE03 

XPNN=XPP 

PR=DE0S00»DE03NN 

IF  (PR.GT.0.D0)  GO  TO  3 

L=L*1 

IF  (L.GT.6)  STOP  1000 

IL  =  0 

XPO=XPOO 

XPN=XPNN 

DEG3O=DE030O 

DE03N=DEO3NN 

XPL=XPO 

XPP=XPN 

DE0SL=DEC30 

DE02R=CEOSN 

XPP=(XPL*XPR)/2.D0 

CALL  STABILUW.VV.I30L) 

V=VV(IKB) 

IF  CDABG(V).GT.l.DO)  STOP  1111 

CALL  ECUILB(V.X,RES.RX,RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RG(6.6.A.WR.W1,MATZ.2.IV1.SV2.IER1) 

IF  (IEP1 .NE.O)  STOP  2222 

CALL  DEG3TBCDE0S.WR) 

DEOSM=DEOS 


TOW00730 
TOW00740 
TOW00750 
TOW00760 
TOW00770 
TOW00780 
TOW00790 
TOW00800 
TOW00810 
TOW00820 
TOW00830 
TOW00840 
TOW00850 
TOW00860 
TOW00870 
TOW00880 
TOWC0890 
TOW00900 
TOW00910 
TOW00920 
TOW009J0 
TOW00940 
TOW00950 
TOW00960 
TOW00970 
TOW00980 
TOW00990 
TOW01000 
TOW01010 
TOW01020 

Towoioro 

TCUOIO^O 
TOW01050 
TOW01060 
TOW01070 
TOW01080 
TOW01090 
TOW01100 
TOW01110 
TOW01120 
TOW01 130 
TOW01140 
TCW01150 
TOW01160 
TOW01170 
TOW01180 
TOW01190 
TOW01200 
TOW01210 
TOW01220 
TOW01230 
TOW012C0 
TOW01250 
TOW01260 
TOW01270 
TOW01280 
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XPM=XPP 

PRL=DE0GL»DE03M 

PRR=DEOSR»DEOSM 

IF  (PRL.GT.O.DO)  GO  TO  5 

XPO=XPL 

XPN=XPM 

DEOSO=DEOSL 

DEOSN=DEOGM 

IL=IL*1 

IF  tlL.GT. ILMAX)  STOP  3100 

DIF=DABS(XPL-XPM) 

IF  (DIF.GT.EP3)  GO  TO  6 

XP=XPM 

GO  TO  C 

IF  (PRR.GT.O.DO)  STOP  3200 

XPO=XPM 


TOW01290 
TOW01300 
TOW01310 
TOW01320 
TOW01330 
TCW013'0 
TOW01350 
TOW01360 
TCW01370 
TOW01380 
TOW01390 
TOWOKOO 
TOW0K10 
TOH01420 
TOW01<;3C 
TOW01440 


XPN=XPR 
DE0S0=DE03M 
DEOCN=DEOCR 
IL=IL-1 

IF  (IL.GT. ILMAX)  STOP  3100 
DIF=DABC(XPM-XPR  ) 
IF  (DIF.GT.EPS)  GO  TO  6 
XP=XPM 
«      LLL=10'L 

NRITE  (LLL.10)  XP.YPP 
3      XPOO=XPNN 

DEOGOO=DEOCNN 
2    CONTINUE 
1  CONTINUE 
STOP 
10  FCRMAT  (2D20.  10) 
1001  FORMAT  ('  ENTER  RANGE  OF  Vp  VARIATION') 

ENTER  NUMEER  OF  INCREMENTS  IN  Yp ' ) 
ENTER  RANGE  OF  Xp  VARIATION') 
ENTER  NUMBER  OF  INCREMENTS  IN  Xp ' ) 
ENTER  EQUILIBRIUM  NUMBER') 
ENTER  DEGREE  OF  STABILITY  CONTROL') 


1002  FORMAT  ( 

1003  FORMAT  ( 

1004  FORMAT  C 

1005  FORMAT  ( 

1006  FORMAT  ( 
2001  FORMAT  (215) 

END 
C 

SUBROUTINE  DSHDEOS.WR) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DIMENSION  WR(6) 

DEOS1=-1.D30 

DO  1  1=1.6 

IF  (WR< I J.LT.DE031 )  GO  TO  1 

DE031=WR( I ) 

IJ=I 

1  CONTINUE 
DEOC2=-1.D30 
DO  2  1=1.6 

IF  (IJ.EQ.Il  GO  TO  2 

IF  (WR(I).LT.DEOS2>  GO  TO  2 

DE0S2=WRC1 ) 

IJJ  =  I 

2  CONTINUE 
DEC3=-1.D30 
DO  3  1=1.6 

IF  (I.EQ.IJ.OR.I.EO.IJJ)  GO  TO  3 
IF  (WR(I).GE.DEOS)  DEOS=WRCI) 

3  CONTINUE 
RETURN 
END 


TOW01450 
TOW01460 
TOW0K70 
TOW01480 
TOW0K90 
TOW0I500 
TOW01510 
TOW01520 
TOKO  15 30 
TOW01540 
TOW01550 
TOW01560 
TOW01570 
TOW01S80 
TOW01590 
TON01600 
TOW01610 
TOW01620 
TOW01630 
TCW01640 
TOW01650 
TOW01660 
TOW01670 
TOW01680 
TOW01690 
TCW01700 
TOW01710 
TOW01720 
TOW01730 
TOW01740 
TOW01750 
TOW01760 
TOW01770 
TOW0178  0 
TOW01790 
TOW01800 
TOW01810 
TOW01820 
TOW01830 
TOWOISSO 
TOW01850 
TON01860 
TOW01870 
TOW01880 
TOW01890 
TOM01900 
TOW01910 
TCW01920 
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FILE:  T0WBIF3   FORTRAN   Al 


PROGRAM  TOWBIF3.FTN 

BIFURCATION  ANALYSIS  OF  TOWING  SYSTEMS 

PARAMETERS  ARE:  Xp.  Lw 

IT  NEEDS  SUBROUTINES  FROM  TOWING. FTN 

USER  DEPENDENT  SUBROUTINES: 

DEGSTB  =  CURVES  ENCLOSING  REGICN  II  OF  FIGURE  13 

(SUBROUTINE  DEGSTB  IS  IN  TOWING. FTN) 
D31     =  CURVES  ENCLOSING  REGION  V  OF  FIGURE  13 

(SUBROUTINE  DS1  IS  IN  SPMBIF) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  MA3SP .NVD. NV . NRD. IZZ . NR . LB .LEN . LW . 

1  NO . NOU . NOUU . NVVV ■ NVRR . NVDD . NVU . NVUU . NRRR . NRVV . 

2  NRDD . NRU . NRUU . ND . NDDD . NDVV , NDRR . NDU . NDUU . NVRD 

DIMENSION  IV1(6).A16.6).VV(3).X(6).WR(6).WI ( 6 ) . Z ( 6 . 6 ) , SVC ( 6 ) 

COMMON/INTGR/ISKEG.NREDP . ITYS. ID- IFDS . ISTAB, IPROP 

CC'MON/SPAR/MASSP .LW.XPP .YPP ,LB 

COMMON/SURGE/SU(7) 

COMMON/XSURG/XU.XUU.XUUU 

COMMON/SWAY/SW( 15 ) 

COMMON/ YAW/YA( 16 ) 

COMMON/MTER/VCAR ,RHO . ABS . CON1 ,CON2 

CCMMON/RE3IST/VEL(<;0  )  .RES I  (40) 

COMMON/VELE/UELU00) 

COMMON/POSTN/X1 -Yl .Zl 

CCMMON/GEOM/AL.RW.G.AET.HW.HW1 

COMMON/PROP/ALE . P . EY . D I  A . ANI U 

COMMON/CTNR/XC(99).YC(99)>2C(99).TC(99) 

COMMON/INT1/IC 

CCMMON/DOC/UC. ALPHA 

CCMMON/UEPT/RLX.RLY.RLZ 

COMMCN/SLOPE/PDRXX.PDRXY.PDRYX.PDRYY 

COMMON/SLAN/RXX6 .RYYb .RXX.RYY 

OPEN    (UNIT^S.FILE^TANKERZ' .STATUS= " OLD' ) 
OPEN    (UNIT=1.FILE=,RE31R' ,STATUS= ' NEW ) 
OPEN    (UNIT=2.FILE=,RES2R1  .  STATUS=  "  NEW  ) 
OPEN    (UNIT=3.FILE=,RES3R> .STATUSi'NEW ) 
OPEN    <UNIT=4.FILE=,RE34R' .STATU3=  "NEW  ) 

CALL  INPUTDC10) 

VCAR  =VCAR»1 . 689D0 

AL    =LW«LB"0.30<;8D0 

ALE   =AL 

MATZ  =0 

IFLOW=l 

1LMAX=1500 

EPS   =l.D-5 


WRITE  (' 

READ  (' 

WRITE  C 

READ  (' 

WRITE  (' 

READ  (• 

WRITE  (• 

READ  (' 

WRITE  (' 

REA3  (• 


1001) 

»)    ai.a: 
1002  ) 

»)     NUMl 
1003) 

n)    bi-e: 

1004) 

« )     NUM2 

1005) 

«)     1KB 


TOW00010 
TCW00020 
TCW00030 
TCW00040 
TOW00050 
TOW00060 
TOW00070 
TOW00080 
TOW00090 
TOW00100 
TOWC0110 
TOW00120 
TOW00130 
TOW00140 
TOW00150 
TOW00160 
TCW00170 
TOW00180 
TOW00190 
TOW00200 
TOW00210 
TOW00220 
TOW00230 
TOW00240 
TCW00250 
TOW00260 
TOW00270 
TCW00280 
TCW00290 
TOW00300 
TCW00310 
TOW00320 

towoo::o 

TOW00340 
TOW00350 
TOW00360 
TCW00370 
TOWC0380 
TOW00390 
TOW00400 
TOW00410 
TOW0C420 
TOW00430 
TOWOO'i'lO 
TOW00C50 
TOW00160 

towoo<;7o 

TOW00480 
TOW00<i90 
TOW00500 
TOW00510 
TOW00520 
TOW00530 
TOW00540 
TOW0C550 
TOW00560 
TOW00570 
TOW005S0 
TOW0C5°0 
TOW00600 
TOW00610 
TOW0C620 
TOW00630 
TOW00640 
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WRITE  ( » ■ 1 006 ] 
READ   («.»)     ID3 
DO  1  1=1. NUM] 

WRITE  O.2001)  I.NUM1 

LW  =  A1*(A2-A1)»(1-1 )/(NUMl-l  ) 

AL  =  LW-LB«0.30<;8D0 

ALE=AL 

XPP=B1 

CALL  STABIL(IVV.VV.ISOL) 

SET  V=VV(K)  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 

V=VV(IKB) 

IF  (DABC(V).GT.l.DO)  STOP  1111 

CALL  EOUILB(V.X.RES.RX.RY) 

CALL  LINEARtX.RES.A.RX.RY) 

CALL  RG(6-6-A.WR.WI.MAT2.r.IVl.SV2,IERl) 

if  (ieri.ne.o)  stop  2222 

if  (ids. ec. 1)  call  degstb ( decs . wr ) 

if  tids.eq.2)  call  dshdeos.wr) 

deosoo=decs 
xpoo  =xpp 

L     =0 

DO  2  J=2.NUM2 

XPP=B1»CB2-B1 )»(J-1)/(NUM2-1) 

CALL  STABILCIVV.VV. ISOL) 

SET  V=VV(K)  FOR  BIFURCATICN  ANALYSIS  OF  K-TH  EQUILIBRIUM 

V=VV(IKB) 

IF  (DABS(V).GT.l.DO)  STOP  1111 

CALL  EOUILB(V.X.RES.RX.RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RGU.6.A.WR.WI  .MATZ.2-IV1  .SV2.IER1  ) 

IF  (IERI.NE.O)  STOP  2222 

IF  (IDS.EO.l)  CALL  DEGSTBtTEOS.WR) 

IF  (IDS.E3.2)  CALL  DSKDEOS.WR) 

DE0SNN=DEC3 

XPNN=XPP 

PR=DEOSOO«CEOSNN 

IF  (PR.GT.0.D0)  GO  TO  3 

L  =  L*1 

IF  (L.GT.6)  STOP  1000 

IL  =  0 

xpo=xpoo 

XPN=XPNN 

DEOCO=CEOSOO 

DECSN=DEOSNN 

XPL=XPO 

XPR=XPN 

DEOSL=DEOSO 

DEOSR=DEOSN 

XPP=(XPL*XPR)/2.D0 

CALL  STABILCIVV.VV. ISOL) 

V=VV(IKE) 

IF  (DABS(V).GT.l.DO)  STOP  1111 

CALL  EOUILB(V-X.RES.RX.RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RGC6.4.A.WR.WI , MAT2 .2 . 1 VI . SV2 . IER 1 ) 

IF  (IERI.NE.O)  STOP  2222 

CALL  DEGSTB(DEOS.WR) 

DEOSM=DECG 

XPM=XPP 


TOW0C650 
TOW00660 
TOW00670 
TOW00680 
TOW00690 
TOW00700 
TOW00710 
TOW00720 

TOW00730 
TOW00740 
TOW00750 
TOW00760 
TOW00770 
TOW00780 
TOW00790 
TOW00800 
TOW00810 
TOW00820 
TOW00830 
TOW008<;0 
TOW00850 
TOW00860 
TOW00870 
TOW00880 
TOW00890 
TOW00900 
TOW00910 
TOW00920 
TOW00930 
TOW00940 
TOW00950 
TOW00960 
TOW00970 
TOW00980 
TOW00990 
TOW01000 
TOW01010 
TOW01020 
TOW01030 
TOW010C0 
TOW01050 
TOW01060 
TCW01070 
TOW01030 
TOW01090 
TOW01100 
TOW01110 
TOW0112C 
TOW01130 
TOW011S0 
TOW01 150 
TOWOlloO 
TOW01170 
TOW01180 
TOW0119O 
TOW01200 
TOW01210 
TOW01220 
TOW01230 
TOWC1260 
TOW0125C 
TOW012tO 
TOW01270 
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PRL=DE03L«DEOGM 

PRR=DE03R«DE0SM 

IF  (PRL.GT.O.DO)  GO  TO  5 

XPO=XPL 

XPN=XPM 

DE0G0=DE03L 

DE03N=DE03M 

IL=IL*1 

IF  (IL.GT.ILMAX)  STCP  3100 

DIF=DAB3(XPL-XPM) 

IF  (DIF.GT.EP3)  GO  TO  6 

XP=XPM 

GO  TO  <, 

IF  (PPR.GT.O.DO)  STOP  3200 

XPO=XPM 

XPN=XPR 

DE030=DE0SM 


TOW01280 
TOW01290 
TOW01300 
TOW01310 
TOW01320 
TOW01330 
TOWOli&O 
TOW01350 
TOW01360 
TOW01370 
TOW01380 
TOW01390 
TOWOUOO 
TOW01610 
TOW01420 
TOW01430 

towoi«<;o 


DEOSN=DEOSR 
IL=IL*1 

IF  ( IL.GT. ILMAX)  STOP  3100 

DIF=DAB3(XPM-XPB) 

IF  (DIF.GT.EP3)  GO  TO  6 

XP=XPM 
i  WRITE  (L.10)  XP.LW 

3     XPOO=XPNN 

DE03OC=DEOSNN 
2    CONTINUE 
1  CONTINUE 
STOP 
10  FORMAT  C2D20. 10 ) 

1001  FORMAT  C  ENTER  RANGE  OF  Lw  VARIATION') 

1002  FORMAT  C  ENTER  NUMEER  OF  INCREMENTS  IN  Lw'} 

1003  FORMAT  C  ENTER  RANGE  OF  Xp  VARIATION') 

1004  FORMAT  C  ENTER  NUMEE3  OF  INCREMENTS  IN  Xp') 

1005  FORMAT  C  ENTER  EQUILIBRIUM  NUMBER') 

1006  FORMAT  ('  ENTER  DEGREE  OF  STABILITY  CONTROL') 
2001  FORMAT  (215) 

END 
C 

SUBROUTINE  D31CDEOS.1-S) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DIMENSION  WR16) 

DE031=-1 .D30 

DO  1  1=1.6 

IF    CWR(I J.LT.DEOS1 )    GO    TO    1 

CE0S1=WR(I ) 

IJ  =  I 

1  CONTINUE 
DEO32=-l.D30 
DO   2    1=1-6 

IF    (IJ.EO.I)    GO   TO  2 

IF    (WR(I ).LT.DEC32)    GO    TO    2 

DEOS2=WRl I ) 

1JJ=I 

2  CONTINUE 
DEOS=-l .DSO 
DO    I    1=1.6 

IF    (I.E0.IJ.0R.I.E3.IJJ)    GO    TO    3 
IF    CWRd  l.GE.DECS)    DEOS  =  WR(I) 

3  CONTINUE 
RETURN 
END 


TOW01450 

TCW01460 

TOU01470 

TCW01480 

TOW014°0 

TOW01500 

TOW01510 

TOW01520 

TOW01530 

TOW01540 

TOW01550 

TOW01560 

TOW01570 

TCW01580 

TOW01590 
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